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ABSTRACT 

Motivated by debris disk studies, we investigate the gravitational microlensing of background 
starlight by a planetesimal disk around a foreground star. We use dynamical survival models to 
0^ ■ construct a plausible example of a planetesimal disk and study its microlensing properties using es- 

' tablished ideas of microlensing by small bodies. When a solar-type source star passes behind a plan- 

\ etesimal disk, the microlensing light curve may exhibit short-term, low-amplitude residuals caused 

by planetesimals several orders of magnitude below Earth mass. The minimum planetesimal mass 
probed depends on the photometric sensitivity and the size of the source star, and is lower when the 
planetesimal lens is located closer to us. Planetesimal lenses may be found more nearby than stel- 
lar lenses because the steepness of the planetesimal mass distribution changes how the microlensing 

■ signal depends on the lens/source distance ratio. Microlensing searches for planetesimals require es- 
sentially continuous monitoring programs that are already feasible and can potentially set constraints 
on models of debris disks, the supposed extrasolar analogues of Kuiper belts. 

^ , Subject headings: gravitational lensing — methods: data analysis — Kuiper belt — planets and 

■ satellites: general — debris disks 

Qh' 1- INTRODUCTION 

Q i The advent of infrared (IR) observatories, such as Spitzer, has revealed a large number of old (> 10^ yr) disks 

■ with IR excesses, known as "debris disks" (Zuckerman 2001; Wyatt 2008). It is generally believed that the excess IR 
] emission emanates from dust produced by collisions between planetesimals in these gas-poor disks. Despite efforts to 

C3 . hnk the IR emission, dust grains and planetesimal populations (Krivov et al. 2008), the properties of the planetesimals 
remain poorly understood. Fundamentally, IR observations only probe the coUisional cascade of dust grains and 
' are mute on the mass distribution of the primordial planetesimals (Heng & Tremaine 2009). It is fair to say that 
^ : no robust constraints have been set on the planetesimal population from studies of debris disks. In fact, how to 
determine the masses of the largest planetesimals remains an open question (Wyatt & Dent 2002). To date, the only 
' empirical constraints on planetesimals come from our own Kuiper Belt (Luu & Jewitt 2002). This provides incentive 
, to develop alternative ways to detect planetesimals, in order to better understand the planetesimals themselves and 

■ also to enhance the scientific impact of planned surveys of debris disks (e.g., the SEEDS survey by Subaru). 

' Gravitational microlensing (e.g., Paczyhski 1996) is a rapidly maturing field that offers such an alternative. It is 
\L} , already an established way of detecting planets as is evident from recent discoveries (Bond et al. 2004; Udalski et al. 
^ 2005; Beauheu et al. 2006; Gould et al. 2006; Bennett et al. 2008; Gaudi et al. 2008; Dong et al. 2009) and theoretical 
^ work (e.g., Mao & Paczyhski 1991; Gould & Loeb 1992; Wambsganss 1997; Griest & Safizadeh 1998; Gaudi et al. 
. . : 1998). Microlensing by small bodies has also been discussed (Bromley 1996; Agol 2002, 2003; di Stefano & Scalzo 
J> 1999a); various possibilities include microlensing by systems with multiple planets (di Stefano & Scalzo 1999b), wide- 
separation planets (Han et al. 2005), extrasolar moons (Bennett & Rhie 2002) and Earth-like moons around ice giants 
(Han 2008). 

. Drawing on and extending some of these established ideas, we suggest that gravitational microlensing provides an 
■ - ' attractive new way to study planetesimal disks that comprise a population of planetesimals with a distribution of masses. 
While the instantaneous probability of microlensing is dominated by the most massive planetesimals in the disk, the 
total number of microlensing events during a sustained observational campaign is more sensitive to intermediate-mass 
planetesimals. Microlensing can probe not only systems like the observed debris disks — a subset of planetesimal disks 
that tend to be large and dynamically hot (Hcng & Tremaine 2009) — but also disks that are dynamically cold and/or 
small. By combining microlensing and IR observations of hot disks, one can potentially extract information about 
total disk masses and planetesimal mass distributions."* With good constraints on planetesimal mass distributions 
at the low- and high-mass ends, one can infer disk ages (Pan & Sari 2005). In general, detecting planetesimal disk 
microlensing events will provide valuable new constraints on models of debris disks. 

The novelty of our approach stems from uniting dynamics and microlensing. Since the formation of planets and 
planetesimal belts from first principles is fraught with uncertainties (Goldreich et al. 2004), we use a simple survival 
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* While absolute masses can be difficult to determine in microlensing (e.g., Gaudi 2002; Han et al. 2005), the mass ratios between bodies 
that contribute to a microlensing light curve can be measured well. 
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model to construct a plausible planetesimal disk in we focus on one example of a large, hot disk. In ^ we tap 
into some ideas from the microlensing literature, use them to discuss the basic phenomenology of microlensing by 
planetesimal disks, and give simple but realistic estimates of the occurrence of disk microlensing events. We argue 
that long-term monitoring for low-magnification events will inform us about the properties of planetesimals that are 
otherwise invisible to conventional methods of detection. The implications of our results and opportunities for future 
work are discussed in 21 

2. DISK PROPERTIES 

2.1. Disk Parameters from a Survival Model 

A planetesimal disk is characterized by its age (tage), mass (Mdisk) and semi-major axis (a), and is composed of 
planetesimals with a range of masses m and radial velocity dispersions ar- To have an age of tagc, a disk must 
have survived all dynamical processes during that time (Heng & Tremaine 2009). Planetesimal disks may be broadly 
separated into those for which planetesimal orbits do or do not cross within their lifetime, respectively termed "hot" 
or "cold" disks. For the purpose of discussion, we adopt the following disk parameters: 

iagc = 10Vr, Mdisk -lOMe, a=10AU. (1) 

Our assumption for a is plausible because dozens of debris disks have been found with 1 < a < 100 AU (see Figure 7 
of Wyatt 2008). The assumed age is also a factor 10 longer than conceivable time scales for protoplanetary disks to 
disperse their gas (Hillenbrand 2008; see also Figure 2 of Wyatt 2008 and references therein). In addition, we assume 
the internal density of the planetesimals to be pp = 3 g cm~'^ and the mass of the parent star to be Mi, = Mq. 
The disk mass can be expressed as 

naaut 

Mdisk = 271 E(a') a' da', (2) 

J ain 

where E is the mass surface density, and the inner and outer disk radii are ain = a/rj and aout — cii]. We adopt 

S(a') « ^-j , (3) 

as is assumed for models of the minimum mass solar nebula (Weidenschilling 1977). Such a scaling also implies that 
the Toomre parameter, Q = ar^/nGYi, is independent of the semi-major axis, where Q = ^/GJdZ[a^ is the orbital 
angular velocity. Equation ([2]) can be rewritten as 

Mdisk = /m7ra2l](a), (4) 

and setting = 1 then implies r/ w 1.28. In other words, our fiducial disk has a radial width of Cout — lin ~ 0.5a, 
and an area Adisk = '"'(aout ~ ^In) ~ t^o?'. 

The maximum allowed planetesimal mass and radial velocity dispersion are mainly determined by the requirements 
that the timescale for gravitational scattering exceeds iagc and that the disk is thin. The scattering condition reads: 



d\nielV 



dt 



> (5) 



where cq is the root mean square eccentricity of the planetesimal orbits. The left-hand side can be evaluated using 
equation (3.29) of Stewart & Ida (2000), where we take the root mean square inclination of the planetesimal orbits to 
be io = eo/2. 

Taking h to be the root mean square disk height, the thin disk condition reads: 

h < 0.35/ea, (6) 

which can be interpreted as the condition that most of the planetesimals are bound or that the radial velocity dispersion 
is less than the circular speed. The preceding expression is equivalent to eo < /e or < we adopt /e = 0.5 

(Heng & Tremaine 2009). 
Equations (O and ^ together yield: 



m < 2Mgj 

ffr < 3.4 km s^' (w 0.4 afl) 



(7) 



For comparison, the circular speed at a = 10 AU is ail « 9.4 km s ^, while the typical bulk velocities of stars near 
the Galactic bulge are ~ 100 km s-\ Thus, as a first approximation we may consider that microlensing is driven by 
the bulk motions of the lens and source, and neglect the motions of the planetesimals within the disk. 

^ For illustration wc consider source stars near the bulge because they are common targets for microlensing campaigns, but as we shall 
see it will bo possible to consider other source locations without an appreciable loss of signal. 
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For our fiducial disk parameters, Toomre stability is trivially fulfilled. Note that for m w 2M(^ and cr^ ~ 0.4 afi, the 
Safronov number is 3> 1. Therefore, i^gc -c^ ^ oc 9^ is a stronger condition than f^^^ > cx 6, where tc is the 
collision time. 

Further details of dynamical survival models for long-lived planetesimal disks are described in Heng & Tremaine 
(2009). 

2.2. Planetesimal Population 

We assume the planetesimals have a mass distribution such that the number of planetesimals between mass m and 
m + dm is 

dN 

— dm = Bm-°' dm. (8) 
dm 

Alternatively, one can instead write the size distribution, dN/dr oc r~*, where g = Sa — 2 and r = (3m/47rpp)^/'^. 
Dohnanyi (1969) showed that a = 11/6 (g = 7/2) in a steady-state system. Pan & Sari (2005) rederived and generalized 
the results of Dohnanyi (1969), allowing for collisions to be inefficient (i.e., the kinetic energy of the bullet does not 
entirely go into breaking up the target); they inferred that 23/8 < q < 22/7 or 13/8 < a < 12/7. While the mass 
function may not be a single power law over the full range of interesting masses (as in our Kuiper Belt; Bernstein et al. 
2004), in this pilot study we use a single power law. If we restrict ourselves to a < 2, and we assume that the smallest 
planetesimal mass in the distribution is much smaller than the largest mass tol, then the normalization factor is 

B« (2-a)Afdiskm^"'. (9) 

For the purpose of microlensing, one needs to consider the projected planetesimal number density per unit mass, 
dn-p-cojl dm, such that 

dN _ ^_ f"-""^ drip, 

from which it follows that 



^ 27r / a' da', (10) 

dm L. dm 



cJjlproj Bm 



dm 



'3/2 

(11) 



3. DISK MICROLENSING 
3.1. Basic Picture 

Figure [T] shows examples of a source star passing behind a fiducial disk. If the source star passes close (in projection) 
to the parent star we get a classic stellar microlensing event, and if it also passes close to a planetesimal we can get a 
short secondary feature in the light curve. This phenomenon is familiar from planetary microlensing (Mao & Paczyhski 
1991; Gould & Loeb 1992; Han et al. 2005); the main difference with planetesimal disk microlensing is that the masses 
are smaller and we must explicitly consider a significant population of planetesimals. 

Qualitatively, there are three limiting regimes of microlensing involving a planetesimal disk: (1) the source star passes 
directly behind a planetesimal that is "far" from its parent star; (2) the source passes directly behind a planetesimal 
that is "near" its parent star; or (3) the source star passes directly behind the parent star in such a way that the light 
curve is still sensitive to the presence of the planetesimal disk. In the first two cases, "far" and "near" are defined 
with respect to the Einstein radius of the parent star. 



i?K. = ^V^ = 4.0AU '''''' 



c 



0.25 8 kpc Mg 



(12) 



where D = DiDis/Dg given the distance to the lens (Di), the distance to the source (Ds), and the distance from the 
lens to the source {Dis = Dg — Di). It is convenient to define the lens/source distance ratio 

- ^, (13) 

and then write D — fi{l~ fi)Ds- The distinction between the "far" and "near" regimes arises because the star creates 
a tidal shear 7 = (i?E*/d)^ at a projected distance d, which enhances the cross section for microlensing by ~ 7^. 
In the "far" regime, we can neglect the influence of the parent star on microlensing by the planetesimal, and the 
system can be regarded as a wide-separation binary (di Stefano & Scalzo 1999a; Han et al. 2005). By contrast, in the 
"near" regime the microlensing signal is significantly affected by "planetesimal caustics" (i.e., analogues of planetary 
caustics; see Schneider & Weiss 1986 for a full discussion of binary point-mass lens systems). Finally, the third regime 
corresponds to a high-magnification event in which the light curve may be perturbed by secondary caustics that can 
be used to detect the presence of planets (Wambsganss 1997; Gaudi et al. 1998; Griest & Safizadeh 1998), or even 
planetesimals. 
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Fig. 1. — The top panel shows a schematic diagram (not to scale) of a source star passing behind a fiducial planetesimal disk with 
-'^disk = 10 Mfg, a = 10 AU, Ds = 8 kpc, /; = 0.1, oil = 2 Af^ and a = 11/6. The thin red circle indicates the Einstein radius of the 
parent star (of mass Mq; -Re* = 1-2 AU), while the thick brown ring represents the planetesimal disk. The dotted lines show sample 
source trajectories that lead to the light curves shown in the middle and bottom panels. In the middle panel, microlensing by the star is 
weak (<5 ~ 1%) compared to that by the planetesimal (<5 ~ 10%). In the bottom panel, the star dominates the microlensing lightcurve, 
while the planetesimal contributes a 5 ~ 1% residual. For illustration, we have assumed v± = 100 km s^^ for the purpose of computing t, 
the time of passage. The insets zoom in on the light curves associated with the planetesimals; tpcak is the time at which the planetesimal 
microlensing events peak. Note that the light curves are from full numerical calculations that cio not involve the assumptions of isolation 
and no shear. While this schematic serves to illustrate the basic idea of planetesimal disk mic rolen sing, the brightness of the lens star for 
nearby events renders, for example, a 10% deviation into a ~ 0.1% deviation for /; =0.1 (see i|4.1l l. 
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The spatial scale for microlensing by a planetesimal is set by its Einstein radius, 

Re - -VCmD « 10^° cm ( ^| . (14) 

c \^0.01Me2 kpc^ ^ ' 

For comparison, a planetesimal of mass m = 0.01 A/0 and density pp = 3 g cm~'^ has a physical size of r '--^ 10* cm, 
implying that it can be treated as a point lens; we return to this issue in Figure 3] Note that since the Einstein radius 
scales as Re oc to^/^ while the physical size scales as r c>c m^/^, there is some mass threshold below which the physical 
size is comparable to or larger than the Einstein radius. When the two are comparable, it may be possible to infer the 
size from combined microlensing and occupation measurements (Bromley 1996; Agol 2002). 
A basic time scale for microlensing is the time it takes to cross the Einstein diameter, 

i^,^^«35min ^V^V (15) 

' vi_ \^0.01Afe2kpcy VlOOkms-V ' ^ ^ 

where v± is the transverse relative velocity of the lens and source. As discussed below, the actual duration of a 
microlensing event may be longer than tE,d because of finite source effects. In particular, as the lens mass decreases 
and iE,d becomes shorter, finite source effects set a floor on the event duration given by the time it takes the lens to 
cross the diameter of the source star, 

2/i^* . I fi R* 



t^d = ^^«116min N^— ^ , (16) 

vi_ \Q.5RqJ VlOOkms-V 

where the factor of /; handles the projection of the source size into the lens plane. The last time scale of interest is 
the time taken to cross the full planetesimal disk, 

tdis.-lyr (l^)(iookm,-i) • (17) 
3.2. Figures of Merit 

While some authors have used detailed simulations of microlensing campaigns to forecast the detection of planets 
(e.g., Han et al. 2005), we elect to use simple but reasonable figures of merit to give initial estimates of the occurrence 
of planetesimal disk microlensing. Wherever possible, we compare our results with those in the literature. We 
ignore factors of order unity related to the planetesimal disk geometry and assume a face-on disk. In this paper, we 
make two assumptions: (1) the planetesimals are far enough from the parent star that we can neglect shear; (2) the 
planetesimals are effectively isolated from one another, so the light curve is influenced by only one planetesimal at a 
time. The assumption of isolation can be verified a posteriori by checking that the microlensing optical depth is low. 
Relaxing these assumptions will provide interesting opportunities for follow-up work, as we discuss in 21 

Microlensing increases the observed flux of the source star by the time-dependent magnification factor A{t) 
(Paczyhski 1986). In general, the total observed flux may include contributions not only from the amplified source star 
but also from the lens star and any other stars along the line of sight that are projected within the same resolution 
element: 

Fohs{t) = A{t)Fsourcc + F\cns + -Pother • (18) 

We set Ficns = -Pothers = in the present analysis for clarity and note that the effects of stellar blending (di Stefano & 
Esin 1995; Han et al. 2006) must be kept in mind when considering prospects for detecting low-amplitude microlensing 
events (see SHI)- In the absence of blending, the maximum fractional change in the fiux is 



6 = max 
t 



Fohs{t) ^ 



max[^(t)]-l. (19) 



For the purpose of gaining a general understanding of planetesimal disk microlensing, we consider a microlensing event 
to be detectable if 6 equals or exceeds some detection threshold Sdct- 

A second factor that may be relevant for detection is the duration of the event. Because of finite cadences and a 
desire to have more than one point on the light curve during a planetesimal disk microlensing event, there may be some 
minimum event duration t^i^ that can realistically be detected, allowing for the definition of two distinct regimes: 

t*,d > tmin ■ "amplitude-limited", 
^*.d < ^min ■ "duratiou-limited" , 

since the shortest possible event duration is given by the source crossing time (equation |16|). The duration-limited 
regime occurs when 

/ ?;_L \ I Ri, 
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We have chosen the nominal duration threshold ^min = 10 min based on cadences that can be achieved today (e.g., 
Dong et al. 2009) and those planned for next-generation microlensing campaigns (Bennett et al. 2009; Gaudi et al. 
2009), but the result can be trivially rescaled. An added complication might be to impose a longer duration threshold 
for lower-amplitude events, perhaps to conserve total signal-to-noise. 

In the amplitude-limited regime, planetesimal disk microlensing events have durations equal to or exceeding ^rnin by 
definition, so we do not need to factor the duration threshold into our calculations. In the duration-limited regime, 
some fraction of events (associated with low- mass planetesimals) will have durations shorter than tmim implying that: 
their chances of detection will be reduced due to finite sampling; and even if they are detected the light curve will not 
be temporally resolved, leading to questions regarding whether they are true detections (see ^4.ip . The exclusion of 
these low-mass events will reduce the total expected number of detected events associated with a given planetesimal 
disk. We do not explicitly consider the duration-limited regime in our calculations, because it depends on the details 
of a particular microlensing survey and on the velocity distribution of source stars. Nevertheless, we shall emphasize 
the circumstances under which our conclusions will be affected by a minimum detectable duration. 

Specifying the condition S > Sdet is equivalent to requiring b > be/,, where b is the impact parameter of the source 
star relative to the lens. As a zeroth order estimate, we expect b^ to be comparable to the Einstein radius R^. It may 
be possible, however, to have a detectable perturbation even when &0 exceeds i?E- To allow for this possibility, we 
write the general form 

60 = 0i?E, (22) 

and call (j) a "boost factor." Such a boost factor has been described by di Stefano & Scalzo (1999a), albeit with no 
name. 

A standard figure of merit is the optical depth, which gives the probability that microlensing is detectable at any 
given instant in time. In the limit that microlensing shear by the parent star is unimportant, each planetesimal has a 
"circle of influence" with area nb'^ (equivalent to the "lensing regions" described by di Stefano & Scalzo 1999a), and 
the optical depth is the fraction of the projected area of the disk (Adisk) that is covered by the circles of influence: 

2W ^ a' da') ^ dm 
dm I Adisk 



7r62 dN 



Aa\^v dm 



dm (23) 



B I 6^ m^"" dm. 



_ AGD 
{acf 

We have neglected a factor of order unity related to disk geometry in equation ([23|) . The limits of integration for the 
mass integral are discussed below. 

A second — and perhaps more interesting — figure of merit is the expected number of microlensing events in one disk 
crossing (equation [17\). Conceptually, this is given by the number of planetesimals in a band of width 2b^ centered 
on the chord traced by the source. If we assume the source crosses the full diameter of the face-on disk, we have: 



N= I [2 I " da' ) 2b^ dm 

dm 



2bd, dN , 

— 3— dm (24) 
ira dm 



dm. 




TTac 



The leading factor of 2 in the first expression enters because the source crosses the planetesimal disk during both 
ingress and egress. The number of events will be reduced relative to this estimate by the fraction of the disk diameter 
actually crossed by the source. (Note that while the event rate depends on the relative velocities of the lens and source, 
the number of events in one disk crossing does not.) 

We can combine Af with the number of disk crossings to estimate the total number of disk microlensing events that 
will occur during a sustained survey. Consider a survey that targets a number N^, of source stars for some period 
^survey, Waiting for lenses pass in front (see §4.2l for other strategies). For simplicity, we assume all the source stars are 
at the same distance Ds- Conceptually, we picture each disk tracing out a band of width 2aout and length wxtsurvcy 
We combine that with the three-dimensional number density of lens stars, ni, to compute the fraction of the sky 
covered by the bands, factor in the number of source stars to obtain the number of disk crossings, and then multiply 
by Af to estimate the total number of planetesimal microlensing events during the survey 

A'survoy = / 2aoutV±tsurvniN^,jV dDi. (25) 
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Of particular interest is the scaling with 

^%^ocn,(/OAA(/0, (26) 

where we highlight the fact that both the number density of lens stars and the number of events per disk crossing 
depend on For comparison, consider the number of classic stellar microlensing events in the same survey, which 
has the same form as equation ([25]) but with Oout b^i, = (I)Re* and J\f ^ 1 (since by definition the number of 
microlensing events when the source crosses a lensing region is unity). Since b^^, oc i?^;^ oc fi{l — fi), we have 



* survey 



dfi 



-^ni{f,)^fi{l-fi). (27) 



Finally, we can quantify the distribution of microlensing event durations. While the Einstein crossing time is given 
by equation (fTS]) . the actual event duration also has a factor of the boost: iicns = 2<j)R-E,/v±. This depends on mass 
(through both i?E and 0), so it is useful to compute the weighted mean: 



{hens} — Tr I *ions dm, 




Note that 



(tlcns) ^ 2t 
idisk TI'A/' 



(28) 



(29) 



The limits of integration for the mass integrals are {frimin, "t-l}- The minimum detectable mass m^in results from 
two conditions. The first is that the Einstein radius of the planetesimal exceeds its physical radius, r > R^, which 
is satisfied for all but the lowest values of m and /; considered (see Figure H|) . The second condition requires the 
planetesimal to be "massive enough" to microlens a source that has some finite projected size R^,. We quantify this 
condition by considering the ratio of the angular source size to the angular Einstein radius, which we denote by 
(see ij3.4[) . Collectively, the two conditions yield the minimum mass 



where 



mm{l/p4) <f. (31) 



The minimum value of effectively quantifies the smallest planetesimal mass that can produce a fractional flux 
change 6 > 5det, given the finite size of the source. To be in the amplitude-limited regime, a third condition on m^in 
that involves ^min is applicable, 

m^.n> ^^^^ . [32) 

As discussed after equation (j21[) . we do not explicitly include this condition in our calculations, but we do consider 
how it will affect our conclusions. 

3.3. Zeroth Order Estimates 

To make initial estimates, we require that the impact parameter be smaller than the Einstein radius (i.e., 0=1), 
which amounts to setting the detection threshold to ddct — 0.34 for a small source. We also require that the Einstein 
radius be larger than the projected size of the source (i.e., e = 1). When (j) is independent of the planetesimal mass, 
equations (O and (^5)1 yield dr/dm oc m^~" oc dM^isk/dm. In other words, the optical depth is directly proportional 
to the total disk mass — but insensitive to how the planetesimal masses are distributed — above TOmin- Since the 
total disk mass is dominated by the high-mass end of the planetesimal mass distribution, we have 



4GD „ 

7 ;2*^disk 

(ac) 



1 - 



^min 



mL 



^^Mdisk. (33) 



(ac) 



If we fix the source and vary the distance to the lens, the distance D — //(I — fi)Ds peaks at /; = 0.5. The peak 
can be shifted slightly becaus e TOm in also depends on /;, but our full calculation indicates that optical depth is still 
maximized near /; « 0.5 (see i j3.5[) . The disk parameters in equation ([T]) yield t s:! 5 x 10~^ for // = 0.5. Thus, even 
if there is a source behind a planetesimal disk, the instantaneous probability of disk microlensing is low. 
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The expected number of microlensing events during one disk crossing is 

-a \ 

(34) 




Here the fact that mmin depends on /; becomes important; it shifts the peak in M to values somewhat lower than 
/; = 0.5, and also determines the slope of the J\f versus /; curve. In particular, rrinun oc /;/(! — /;) while l) cx /;(! — /;), 
so for small /; we find 

M cx /f-" = f^/^ (35) 

for a = 11/6. With = 8 kpc, i?* = Rq and ttil = 2Mq, we find that TV peaks when /; w 0.09 (mmin ~ O.OIM®), 
with a value of maxjA/"} « 4 x 10^"^. This is still in the fiducial amplitude- limited regime (cf. equation [2T|) . so the 
scaling with /; is accurate. The scaling differs in the duration-limited regime, but this is of less interest to us in the 
present analysis. 
Using equation ()29|) . the mean microlensing duration is 




(ilcns) « — ( ] ./.-^^ 3/2-0 - (36) 

Adopting the same parameter values as in the preceding paragraph, we have (tions) ~ 2 hours when /; « 0.09. 

3.4. Finite Source Effects 

The simple estimates we just obtained should underestimate the occurrence of planetesimal disk microlensing for two 
reasons. Firstly, the assumed detection threshold of S^ct = 0.34 was very generous. Secondly, we explicitly neglected 
planetesimals whose Einstein radius is smaller than the projected size of the source star. If we account for finite source 
microlensing, we can consider the possibility that smaller planetesimals may contribute to the overall signal (di Stefano 
& Scalzo 1999a). 

In this analysis, it is useful to introduce angular variables: 

6k = 77-, Ob = — , 9^ = — ^, (37) 

i-'l i-'s 

where h is the impact parameter. The source size is important in relation to the Einstein radius, so we consider the 
ratio 

P. ^ 1^. (38) 

If we define 

then the lensing magnification A{uq^ p^,) formally involves a two-dimensional integral, but Heyrovsky & Loeb (1997) 
and Lee et al. (2009) give a useful algebraic approximation (see Appendix), allowing computations to be significantly 
accelerated. 

We evaluate A in the plane of Oe/Os and Ob/Oi,. Some portion of this phase space satisifies 5 > S^ct', this region is 
bounded by the curve, 

i; = i,{mji,b^b^), (40) 

which can be inverted to solve for (p and e. Examples of t/j curves for Sdct = 0.01 and 0.1 are shown in Figure [21 For 
comparison, the shaded region in the figure shows the region of phase space used for the zeroth order estimates in ^3.3\ 
(specifically, the region bounded by the curves 6i, = 9e and b = Re)- Allowing a more sensitive detection threshold 
and accounting for finite source microlensing clearly increases the region of phase space that can contribute to the 
microlensing signal. 

To illustrate the enhancement, we plot the boost factor as a function of the planetesimal mass in Figure The 
shape of the (/>(m) curve is not very sensitive to the value of which can be understood heuristically — the maximum 
impact parameter roughly corresponds to a situation in which the near edge of the source star is some distance CRe 
away from the lens, where C is a constant that depends on the detection threshold (5dct- 

b^^fiR^+CRE, (41) 

where /^i?^ is the size of the source star projected into the lens plane. This is equivalent to the boost factor being 

(42) 

Re 

® This figure is reminiscent of the bottom panel of Figure 6 of Han et ah (2005), where they adopted /; = 0.75 with Ds = 8 kpc. 
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Fig. 2. — In the plane of impact parameter and Einstein radius (both in angular units, normalized by th e siz e of the source star), the 
shaded region denotes the phase space used for zeroth order estimates of the microlensing figures of merit f i|3.3|l . The -0 curves show how 
the allowed phase space is enlarged when we consider different detection thresholds ((5dot) a-nd account for finite source effects ( i|3.4l l. 




Fig. 3. — Boost factor </> as a function of planetesimal mass m, for different values of the lens/source distance ratio (/;) and the detection 
threshold (iSdct)- The dotted vertical lines indicate the minimum detectable mass if the Einstein ra dius is required to be larger than the 
projected source size (e = I). The dot-dot-dot-dash curve shows the heuristic ij> curve from equation l|42|l . 

The constant C is described by equation (4) of di Stefano & Scalzo (1999a); we checked that our values of C in Figure 
[3] are consistent with their formula. Since Re oc m^/^, we expect that is approximately constant at large m, then 
rises as m decreases, down to a minimum threshold mass. This heuristic shape is shown in Figure [31 We see that this 
general argument explains the overall shape of the (j) curves, although it is too simple to capture the full complexity 
near the minimum mass. 

Finite source effects dictate that the boost factor can be as high as ^3 or 12, for a detection threshold of Sdct = 0.1 or 
0.01, respectively (i.e., a ten or one percent perturbation). A second result is that the minimum mass that contributes 
to microlensing can be smaller than one naively estimates by requiring the Einstein radius to be larger than the 
projected source size. For Sdet = 0.1 or 0.01, mmin is reduced (relative to the naive estimate) by a factor of 1/e ~ 20 or 
200, respectively. For example, if the disk is located at /; = 0.5 and we consider ^dct = 0.1, we have mmin ~ 5 x lO^'^Af® 
instead of about 0.1 M^. More generally. Figure S] shows function of /; for different detection thresholds, 

demonstrating that disk microlensing can, in principle, probe planetesimals down to several orders of magnitude below 
Earth mass. 

It is important to note that Han et al. (2005) conclude it will be difficult to probe masses below ^ 0.02 (at 
fi = 0.75) with microlensing. One key difference is that Han et al. (2005) consider a stringent limit for their signal- 
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Fig. 4. — Minimum detectable mass rrimin as a function of the lens/source distanc e ra tio for different values of the detection threshold 
(5det- The shaded region shows the nominal duration-limited regime (see equation |21)'l. The thin dot-dash line shows the point at which 
the Einstein radius is the same as the size of planetesimal (i?E = r); for the detection thresholds we consider, the point lens approximation 
breaks down only at /; < 10"'^. We note that when i?E ~ "r there can be simultaneous microlensing and occultation, which probe the 
physical size of the lens (Bromley 1996; Agol 2002). The horizontal dashed and dotted lines indicate the masses of the Moon and Ceres, 
respectively. 
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Fig. 5. — Figures of merit for planetesimal microlensing as a function of the lens/source distance ratio optical depth t (left); and 
expected number of events per disk crossing A/" (right). The various curves are for different values of t he largest planetesimal mass ttil 
and the detection threshold <5dot- (The black curves correspond to the zeroth order estimates in ^13.31) The curves are reliable in the 
amplitude-limited regime, but are overestimates in the duration-limited regime (shaded; see equation \21\ ). 

to-noise ratio of S/N = VTOOO. For comparison, with /; = 0.75 and ^dot = 0.1 we find mmin 0.02 Mq (Figure [S]). 
Using equation (17) of Han et al. (2005), we obtain S/N ~ 29 for our value of TOmin, assuming only one detection 
point per crossing. Plugging in our estimate for S/N into their equation (18) then yields a minimum mass that is 
consistent with our computed value. A less stringent criterion for the S/N threshold (e.g., S/N w 8 by Gaudi et al. 
2002) will lead to a different value for TOmin- If we consider ^det = 0.01, mmin ~ 2 x 10^'^ (Figure[3]) and retain all 
of the numbers in the preceding estimate, we obtain S/N ~ 6 instead and are consistent with equation (18) of Han et 
al. (2005) to within a factor of 2. Therefore, we believe our calculations are consistent and the issue is more one of 
deciding what constitutes a detectable event (see N4.1l for more discussion). 

3.5. Improved Estimates of the Figures of Merit 

We can now obtain better estimates of the optical depth and number of events per planetesimal disk crossing, taking 
into account finite source effects to properly treat the full range of allowed lens masses. Figure [5] shows examples of 
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Fig. 6. — Mean microlensing event duration (cf. equation 1281 ). as a function of the lens/source distance ratio /;. The horizontal dashed 
line indicates (ticns) = 10 min, which is about the best cadence that current microlensing surveys are capable of. The curves are reliable 
in the amplitude-limited regime, but are underestimates in the duration-limited regime (since events that are too short to be detected are 
excluded). We have assumed v± = 100 km s^^, but the duration can be trivially rescaled for other choices. 

T and A/" as a function of // for different detection thresholds as well as two values of the largest planetesimal mass, 
mL = 0.1 and 2 M^; the latter value is the maximum allowed mass for our planetesimal disk configuration (see §2.ip . 
For illustration, we continue to set a = 11/6. 

The optical depth peaks when the lens is approximately halfway to the source (/; ~ 0.5), and scales as r cx /; when 
the lens is near the observer (i.e., /; is small). Reducing mL from 2 to 0.1 means the disk mass is distributed 
among a larger number of intermediate-mass planetesimals. This has little effect on the zeroth order estimate for the 
optical depth (the black curve in the left panel of Figure [Sj cf. equation [33]), but it does cause a modest increase 
in the more realistic optical depth estimates — at least for small // values — because the boost factor enhances the 
effects of intermediate-mass planetesimals. For Sdct = 0.1, the peak and drop-off in the r curve occur at smaller values 
of /; for mL = 0.1 M0 than for mL = 2 M®, simply because the planetesimals become too small to produce S > Sdct- 
On average, the more realistic optical depth is boosted (relative to the zeroth order estimate) by a factor ~3 or 20 for 
a detection threshold of Sdct =0.1 or 0.01, respectively. Even so, the optical depth remains small. 

The number of events per disk crossing is considerably more sensitive to the value of mL, because this figure of merit 
is more sensitive to low-mass planetesimals. For mL = 2 M^, the realistic estimate for A/" is a factor ~9 or 60 larger 
than the zeroth order estimate, for Sdet = 0.1 or 0.01, respectively; and for mL = 0.1 the enhancement factors 
are ~19 or 140. Nevertheless, we still find Af < 1 for typical parameters, which means we are unlikely to see multiple 
microlensing events from a single disk (for our fiducial disk parameters). 

The mean event duration decreases with /; , as shown in Figure [51 This actually reconciles the low optical depth 
with the modest total number of events: any single observation is unlikely to catch a short, transient event, but a 
thorough monitoring program will be able to find it. 

The weak dependence of Af on /; stems from having a steep mass distribution of planetesimals. If all of the 

1/2 

planetesimals had the same mass, we obtain the scaling Af oc (for small fi). Instead, with a planetesimal mass 
distribution dN/dni oc m^" we find N oc ff " (again for small /;). The dependence on a is illustrated in Figure [71 
which shows Af versus /; for the other a values described in i j2.2l We see that for a < 11/6, Af decreases by factors 
~ 2 to 3, but still remains non-negligible. All of the scalings discussed here apply to the amplitude-limited regime; the 
scalings with /; are steeper in the duration-limited regime, because at low fi the events created by the lowest mass 
planetesimals become too short to be reliably detectable. 



3.6. Distribution of Distance Ratios for Detected Events 

Having understood how the figures of merit depend on /;, we can now make simple estimates of the distribution of 
/; values that may be found in a microlensing survey. The key difference from the preceding analysis is that we need 
to combine the raw figures of merit with a reasonable estimate for the spatial distribution of lens stars. 

We use equations ([26]) and ([27| to interpret the probability density of detecting microlensing events at a given /; 
value as 

/ r \ ^A^survcv 
Pdisk (Jl) oc — , 

(43) 

a\ dN survey,*- 
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Fig. 7. — Expected number of microlensing events per planetesimal disk crossing for various values of the planetesimal mass distribution 
index a. We adopt rriL = 0.1 Mjg and 5dct = 0.01 for illustration. The shaded region again shows the nominal duration-limited regime. 

The proportionality constant in each case serves to normahze the probabiHty density distribution such that it integrates 
(over /z) to unity. 

The source stars are assumed to be in the Galactic bulge. For the lens stars residing in the disk of the Galaxy, the 
density distribution can be approximated with a simple exponential disk model, 

ni(xexp(^-^y (44) 

where Rg — {1 — fi)Ds is the distance from the Galactic center and H = 3 kpc is our adopted value for the radial 
scale length (sometimes termed "radial scale height") of the Galactic disk (Ojha 2001). For the lens stars residing in 
the Galactic bulge, the density distribution is 

m cc (45) 

where the power law index is either s = 1.75 (Binney et al. 1991) or s = 1.85 (Kent 1992). 

Our estimates for pdisk and p^, are shown in Figure[8l where we assume tol = 0.1 and a = 11/6. For bulge-disk 
microlensing, the classic events peak at /; « 0.85, while the planetesimal disk events peak at /; « 0.69 ((5dct — 0.01) 
and // Ri 0.62 (Jdot — 0.1). It is interesting to see that the /; distribution for a survey with the poorer detection 
threshold {Srct = 0.1) is shifted to lower values compared to a survey with the better detection threshold {Srct — 0.01). 
The reason is that for Sdct = 0.1 the minimum detectable mass can exceed the maximum planetesimal mass, i.e., 
"m-min > n^hj for high values of /;, which causes the microlensing signal to vanish about /; « 0.95. For bulge-bulge 
microlensing, all of the events peak at /; > 0.87 for 1.75 < s < 1.85; we adopt s = 1.8 in Figure [8] for illustration. In 
all cases, planetesimal disk microlensing events always peak at lower /; values compared to classic microlensing events. 

We note that Kiraga & Paczyhski (1994) estimate bulge-bulge microlensing events to be about 5 times more likely 
than bulge-disk ones. 

4. DISCUSSION 
4.1. Detectability 

In this paper, we have adopted the simple criterion that a planetesimal disk microlensing event occurs when a 
planetesimal amplifies background starlight above some specified threshold for longer than some minimum duration. 
While this criterion does not account for all of the subtleties involved in a realistic observational campaign, we believe 
it is useful for elucidating the general properties of planetesimal disk microlensing. Conceptually, our main conclusions 
are: 

1. It is conceivable to detect planetesimals in the range of the mass of the Moon or even Ceres — and possibly 
below — although pushing to lower masses becomes increasingly challenging. 

2. Planetesimal disk microlensing events are likely to be shifted to lower values of /; than classic stellar microlensing 
events, because decreasing /; lowers the minimum planetesimal mass that contributes to the microlensing signal, 
thereby increasing the number of available deflectors. 

Quantitatively, we estimate that if it is possible to detect microlensing events in which a solar-type source star is 
amplified by at least 10%, it will only take a few tens of disk crossings to discover one planetesimal microlensing event; 



Planetesimal Disk Microlensing 



13 




Fig. 8. — Probability density for the lens/source distance ratio /; for both classic (p*) and planetesimal disk microlensing (pdisk) events. 
The source stars are assumed to be in the Galactic bulge, while we consider the lens stars to be both in the Galactic bulge (thin lines) and 
disk (thick lines). The nominal duration-limited regime is again shaded. 

and if it is possible to detect a 1% brightening, it will only take a handful of planetesimal disk crossings (for our fiducial 
disk parameters). 

It is still important to consider some practical aspects of microlensing, especially with regards to understanding how 
they will affect our conclusions. By our adopted criterion, a single deviant data point in a light curve is regarded as a 
detection. A cautionary example comes from the study of Sahu et al. (2001), who detected 6 temporally unresolved 
spikes in their data while scrutinizing the globular cluster M22, and initially interpreted them as microlensing events 
caused by free-floating planets with masses ~ 0.25Af j, where Mj is the mass of Jupiter. Gaudi (2002) argued on 
dynamical grounds that it is implausible for free-floating planets to exist in the core of M22, and that any such planets 
in the halo will not be numerous enough to explain the observed microlensing optical depth. Eventually, Sahu et al. 
(2002) discovered that the "microlensing spikes" were caused by cosmic rays. 

The statistical significance of a microlensing event ultimately depends not just on the maximum fractional flux 
deviation but on more sophisticated quantities such as the total signal-to-noise ratio {S/N; Han et al. 2005) or the 
difference between the goodness of fit for model light curves with and without planets/planetesimals (Gaudi & Sackett 
2000). There is no solid consensus on what should constitute a detection and various thresholds have been adopted 
in the literature: S/N = VlOOO w 32 (Han et al. 2005), S/N = y/OO « 8 (Gaudi et al. 2002) or S/N = 12.5 (Bennett 
et al. 2004), and Ax^ > 100 or 225 (Gaudi & Sackett 2000). For our example in gS31 the detection thresholds of 
Sdct = 0.01 and 0.1 are equivalent to S/N ^ 6 and 29, respectively. 

If the planetesimals have a steep mass distribution function, then the number of events per disk crossing remains 
non-negligible even for nearby planetesimal disks (// <C 1). In principle, the closer the disk, the farther down the 
planetesimal mass function we can probe. However, events with small fi are expected to be relatively rare given 
realistic lens and source populations (see Figure [8]). Also, their short durations will make them more difficult to 
detect, since at fixed cadence an event due to a nearby planetesimal event will have fewer deviant light curve points 
and thus a lower total S/N than a corresponding event due to a more distant planetesimal disk. Finally, for events 
with fi<S^l the microlensing signal will be more diluted by light from the lens star. Consider our illustration in Figure 
[U where /; = 0.1. Even without taking dust extinction into account, the lens star is ^100 times brighter than the 
source star. Since the lens and source cannot be resolved during the microlensing event, a deviation of 10% in the 
source fiux will become merely a ~0.1% deviation in the total observed flux. This issue becomes yet more severe for 
systems with smaller /; values. These issues seem to make it very challenging to detect microlensing events from very 
nearby planetesimal disks. Nevertheless, we recall that planetesimal disk microlensing events are still expected to have 
somewhat lower values of /; than stellar microlensing events (Figure [5]) . 

4.2. Observational Strategies 

Traditional microlensing surveys usually target a field of source stars and wait for lenses to pass in front. Next- 
generation campaigns will use wide-field cameras mounted on l-2m telescopes to monitor ~ 10 deg^ fields at cadences 
of 10 minutes (Bennett et al. 2009; Gaudi et al. 2009). Such cadences should be able to detect microlensing events 
associated with all but the closest planetesimal disks. At this point it is difficult to forecast in detail the number 
of planetesimal disk microlensing events that will be detected due to uncertainties in the abundance of planetesimal 
disks. 

In the wide-separation limit considered, a planetesimal disk microlensing event may or may not be accompanied by 
an event associated with the parent star — isolated, short-duration events may occur (as has been discussed for planets; 
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Di Stefano & Scalzo 1999a,b) and will be quite interesting. For an isolated planetesimal event, it may be difficult to 
determine the actual planetesimal mass because the Einstein radius crossing time is a degenerate combination of the 
lens mass, the lens and source distances, and the transverse relative velocity (Gaudi 2002; Han et al. 2005)7 For an 
event that shows effects from both the planetesimal and the parent star, the light curve alone will strongly constrain 
the mass ratio between the two lenses. Working with mass ratios will already enable interesting planetesimal disk 
science. 

A plausible, alternate strategy will be to target known debris disks and wait for source stars to move behind them. 
As we have seen, the microlensing signal is appreciable even for disks that are relatively nearby. The practical challenge 
is to locate a disk with a suitable star behind it. Typical debris disk searches tend to avoid crowded stellar fields, 
because of the difficulties with the point spread function subtraction that impede the measurement of an IR excess 
(A. Moro-Martm 2009, private communication). In this regard it is unclear whether existing debris disk samples offer 
good candidates for microlensing follow-up. It is worthwhile to consider whether there are observational strategies 
that can combine debris disk observations with microlensing to reap the benefits of both. 

While discovering planetesimal disk microlensing events will obviously be exciting, even the non-detection of plan- 
etesimals in microlensing lightcurves will set interesting upper limits on their masses that will be useful to models of 
debris disks. Wyatt & Dent (2002) proposed that dust clumps embedded in the debris disk of Fomalhaut are the result 
of collisions between planetesimals that may be as large as ~ 1000 km in size (~ 0.2 lunar masses or ~ 13 times the 
mass of Ceres) . Such estimates hinge on uncertain extrapolations based on 450 /^m and 850 /im observations of ~ 7 /^m 
and ~ 0.2 m objects. As mentioned in SJH determining the size of the largest planetesimal remains an open question 
in debris disk studies. One can begin to address this question by examining the microlensing lightcurves of hundreds 
of source (dwarf) stars and searching for statistically significant residuals. In principle, if the distances to the source 
and parent stars, the relative velocity between them, and the planetesimal disk geometry are known, one can infer the 
maximum planetesimal mass detectable for a given magnification threshold. Non-detections will also set constraints 
on the size of the planetesimal disk. Surveying a large number of stars for planetesimal microlensing events will shed 
light on the frequency of planetesimal disks with or without planets, thus constraining models of planet formation. 

4.3. Future Work 

In this paper, we have considered large disks in which the planetesimals lie "far" from the parent star (relative to 
the stellar Einstein radius), and focused on events in which the source passes close to a planetesimal. There are two 
interesting ways to extend our analysis. One way is to account for tidal shear, which may be created not only by 
the parent star but also by other planetesimals; this will allow us to analyze small disks. The second approach is to 
consider a scenario in which the source passes so close to the parent star that the light curve is affected by caustics 
created by the planetesimals. This is a direct analog of high-magnification microlensing events that is used to detect 
planets (Wambsganss 1997; Griest & Safizadeh 1998; Gaudi et al. 1998; Gould 2008), but generalized from the case of 
one or a few massive planets to many planetesimals. The caustics are sensitive to the full population of planetesimals, 
so a high-magnification microlensing event will probe the entire planetesimal disk simultaneously. 

A shortcoming of our idealized analysis is the exclusion of blending effects, in particular the assumption in equation 
(I18p that -Fions = 0. In practice, the lens and source stars are located within each other's Einstein radius during 
a microlensing event, and within the resolution limit of current ground- and space-based telescopes. We intend to 
explore this issue in a future paper. 

Clearly there is much fertile ground for further work on both the formal and practical aspects of planetesimal disk 
microlensing. We believe the possibility of obtaining a new way to analyze planetesimal disks will make such studies 
interesting and exciting. 
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APPENDIX 

FUNCTIONAL FORM OF MAGNIFICATION WITH FINITE SOURCE 

The magnification of a finite source by a point lens has been studied by Heyrovsky & Loeb (1997) and Lee et al. 
(2009). Equation (7) from Lee et al. (2009) gives a useful approximation for the case of a uniform, circular source: 



yl(uo,p*) 



, uq > P*, 



(Al) 



^ When ~ the resulting magnification pattern has structure on the scale of the projected source size, which allows to be 
measured, thus breaking some of the degeneracy. 
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where we have defined 



6crit = arcsin (p^/uo), 



(mo + P*) \J [uq + p*f + 4 - (mo - p*) \J (uo - p*)^ + 4 



(A2) 



(9) = M2 Y + 4 - ui ^ 4. 
The quantities ui and U2 are given by 



Ul 



uq cos 6 — y pJ — Wq sin , uq > Pi, and 6 < 6crit, 
, otherwise, 



and 



U2 



uq cos & + y pI ~ Uq sin^ 9 , uq < p* or {wo > C and 9 < 9ciit} : 
, otherwise. 



(A3) 



(A4) 



The accuracy of the formula increases with the resolution, rij-cs- 

More generally, Heyrovsky & Loeb (1997) give formulae for the microlensing of elliptical sources with non-uniform 
surface brightness profiles. 
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